Matrix product state comparison of the numerical renormalization group and the 
variational formulation of the density matrix renormalization group 
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Wilson's numerical renormalization group (NRG) method for solving quantum impurity models 
yields a set of energy eigenstates that have the form of matrix product states (MPS). White's 
density matrix renormalization group (DMRG) for treating quantum lattice problems can likewise 
be reformulated in terms of MPS. Thus, the latter constitute a common algebraic structure for both 
approaches. We exploit this fact to compare the NRG approach for the single-impurity Anderson 
model to a variational matrix product state approach (VMPS), equivalent to single-site DMRG. 
For the latter, we use an "unfolded" Wilson chain, which brings about a significant reduction in 
numerical costs compared to those of NRG. We show that all NRG eigenstates (kept and discarded) 
can be reproduced using VMPS, and compare the difference in truncation criteria, sharp vs. smooth 
in energy space, of the two approaches. Finally, we demonstrate that NRG results can be improved 
upon systematically by performing a variational optimization in the space of variational matrix 
product states, using the states produced by NRG as input. 

PACS numbers: 78.20.Bh, 02.70.+C, 72.15.Qm, 75.20.Hr 



I. INTRODUCTION 

Wilson's numerical renormalization group (NRG) is a 
highly successful method for solving quantum impurity 
models which allows the non-perturbative calculation of 
static and dynamic properties for a variety of impurity 
models.— >^>^>^>^i^ NRG is formulated on a "Wilson chain", 
i.e. a tight-binding fermionic quantum chain with hop- 
ping matrix elements that decrease exponentially along 
the chain as A^"/^, where A > 1 is a discretization pa- 
rameter defined below and n > is the chain's site index. 
It is thus not applicable to real space quantum lattice 
problems featuring constant hopping matrix elements. 
For these. White's density matrix renormalization group 
(DMRG) is the method of the choice.^^"^ It has been 
known for some tim o^°i^^ that the approximate ground 
states produced by DMRG have the form of matrix prod- 
uct states (MPS) (see Eq. ^ below) that had previ- 
ously arisen in certain stochastic models!^ and quan- 
tum information processing.^ This fact can be exploited 
to reinterpret the DMRG algorithm (more precisely, its 
one-site finite-size version) as a variational optimization 
scheme, in which the ground state energy is minimized 
in the space of all matrix product states with specified 
dimensions To emphasize this fact, we shall refer to 
DMRG as "variational matrix product state" (VMPS) 
approach throughout this paper. 

Quite recently it was understood^^ that NRG, too, in 
a natural way produces matrix product states. In other 
words, when applied to the same Wilson chain, NRG and 
VMPS produce approximate ground states of essentially 
the same MPS structure. The two approximate ground 
states are not identical, though, since the two methods 
use different truncation schemes to keep the size of the 
matrices involved manageable even for very long Wil- 



son chains: NRG truncation rehes on energy scale sepa- 
ration, which amounts to discarding the highest-energy 
eigenstates of a sequence of effective Hamiltonians, say 
Tin, describing Wilson chains of increasing length n and 
yielding spectral information associated with the energy 
scale A~"/^. This truncation procedure relies on the ex- 
ponential decrease of hopping matrix elements along the 
Wilson chain, which ensures that adding a new site to 
the Wilson chain perturbs it only weakly. In contrast, 
VMPS truncation relies on singular value decomposition 
of the matrices consituting the MPS, which amounts to 
discarding the lowest- weight eigenstates of a sequence of 
reduced density matrices.— This procedure makes no spe- 
cial demands on the hopping matrix elements, and indeed 
works also if they are all equal, as is the case of standard 
quantum chain models for which DMRG was designed. 

The fact that a Wilson chain model can be treated by 
two related but inequivalent methods immediately raises 
an interesting and fundamental methodological question: 
How do the two methods compare? More precisely, to 
what extent and under which circumstances do their re- 
sults agree or disagree? How do the differences in trunca- 
tion schemes manifest themselves? VMPS, being a vari- 
ational method operating in the same space of states as 
NRG, will yield a lower-energy ground state than NRG. 
However, it variationally targets only the ground state 
for the full Wilson chain, of length TV, say. In contrast, 
NRG produces a set of eigenenergies {EJj} and eigen- 
states for each of the sequence of effective Hamil- 
tonians Hn, with n < N, mentioned above. Prom these, 
a wealth of information about the RG fiow, fixed points, 
relevant and irrelevant operators, their scaling dimen- 
sions, as well as static and dynamic physical properties 
can be extracted. Are these accessible to VMPS, too? 

The goal of this paper is to explore such questions. We 
shall exploit the common matrix product state structure 
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of the NRG and VMPS approaches to perform a sys- 
tematic comparison of these two methods, as apphed to 
the single-impurity Anderson model. It should be em- 
phasized that our purpose is not to advocate using one 
method instead of the other. Instead, we hope to arrive 
at a balanced assessment of the respective strengths and 
weaknesses of each method. 

In a nutshell, the main conclusion (which confirms and 
extends the results of Ref. [ish is the following: when ap- 
plied to a Wilson chain with exponentially decreasing 
hopping, the VMPS approach is able to fully reproduce 
all information obtainable from NRG, despite being vari- 
ationally optimized with respect to the ground state only. 
The reason is that the VMPS ground state is character- 
ized by products of matrices of the form n!^=o ^^'^"^ (de- 
tails will be explained below), where the matrices with 
the same index n contain information about the energy 
scale A~"/^. As will be shown below, this information 
can be used to construct eigenenergies {E^} and eigen- 
states {l-E^)} for a sequence of effective Hamiltonians Hn 
in complete analogy with (but not identical to) those of 
NRG. The agreement between NRG and VMPS results 
for these eigenenergies and eigenstates is excellent quan- 
titatively, provided sufficient memory resources are used 
for both (and A is not too close to 1, see below). In this 
sense, NRG and VMPS can be viewed as yielding essen- 
tially equivalent results when applied to Wilson chains 
amenable to NRG treatment. In particular, all physical 
properties obtainable from the eigenspectra and eigen- 
states of NRG can likewise be obtained from those of 
VMPS. 

Nevertheless, NRG and VMPS do differ in perfor- 
mance, flexibility and numerical cost. Firstly, since NRG 
truncation relies on energy scale separation, it works well 
only if the discretization parameter A is not too close to 
1 (although the continuum limit of the model is recov- 
ered only in the limit A ^ 1). This restriction does not 
apply to VMPS. Indeed, we shall find that NRG and 
VMPS agree well for A = 2.5, but less well for A = 1.5. 
This in itself is not surprising. However it does illus- 
trate the power of VMPS to get by without energy scale 
separation. This very useful feature can be exploited, 
for example, to obtain resolve sharp spectral features at 
high energies in dynamical correlation functions,— using 
projection operator techniques. However, the latter re- 
sults go beyond the scope of the present paper and will 
be published separately. 

Secondly, since VMPS does not rely on energy scale 
separation, it does not need to treat all terms in the 
Hamiltonian characterized by the same scale A~"/^ at the 
same time, as is required for NRG. This allows VMPS to 
achieve a significant reduction in memory cost compared 
to NRG for representing the ground state. To be specific: 
For NRG, we use the standard "folded" representation of 
the Wilson chain, in which each site represents both spin 
down and spin up electrons, with the impurity site at 
one end (see Fig. [l](a) below). However, it turns out 



that apart from the first few sites of the folded chain, the 
spin-down and -up degrees of freedom of each site are ef- 
fectively not entangled with each other at all (see Fig. [3] 
below). For VMPS, we exploit this fact by using an "un- 
folded" representation of the Wilson chain insteadp^ii^ 
in which the spin up and spin down sites lie on opposite 
sides of the impurity site, which sits at the center of the 
chain (see Fig. [UJb) below) . This unfolded representation 
greatly reduces the memory cost, as characterized by the 
dimensions, D for NRG or D' for VMPS, of the effective 
Hilbert spaces needed to capture the low energy proper- 
ties with the same precision: We find that with the choice 
D' = 2"^Vd, VMPS can reproduce the results of NRG 
in the following manner: (i) if m = 0, the NRG ground 
state is reproduced qualitatively; (ii) if m = 1, all the 
"kept" states of NRG are reproduced quantitatively; and 
(iii) if m = 2 all the "kept" and "discarded" states of NRG 
are reproduced quantitatively. However, in cases (ii) and 
(iii) the reduction in memory costs of VMPS is some- 
what offset by the fact that the calculation of the excited 
eigenstates needed for the sake of direct comparison with 
NRG requires diagonalizing matrices of effective dimen- 
sion D'^. Note, nevertheless, that all information needed 
for this comparison is already fully contained within the 
VMPS ground state characterized by dimension D' , since 
its constituent matrices contain information from all en- 
ergy scales represented by the Wilson chain. 

The paper is organized as follows: Section |ll] sets the 
scene by introducing a folded and an unfolded version of 
the Wilson chain. In Sections Hill and HV] we review the 
NRG and VMPS approaches for finding the ground state 
of a folded or unfolded Wilson chain, respectively, em- 
phasizing their common matrix product state structure. 
We also explain how an unfolded MPS states may be 
"refolded", allowing it to be compared directly to folded 
NRG states. In Section [V] we compare the results of 
NRG and VMPS, for ground state energies and over- 
laps (Section [VAp . excited state eigenenergies and den- 
sity of states (Section IVBp . and the corresponding en- 
ergy eigenstates themselves (Section IV Cp . This allows 
us, in particular, to obtain very vivid insights into the 
differences in the truncation criteria used by the NRG 
and VMPS approaches, being sharp or smooth in en- 
ergy space, respectively (Figs.lSlto fTOl) . In Section [Vll we 
demonstrate that NRG results for the ground state can 
be improved upon systematically by first producing an 
unfolded "clone" of a given NRG ground state, and sub- 
sequently lowering its energy by performing variational 
energy minimization sweeps in the space of variational 
matrix product states. Finally, Section [VTTI contains our 
conclusions and an assessment of the relative pros and 
cons of NRG and VMPS in relation to each other. 
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II. FOLDED AND UNFOLDED 
REPRESENTATIONS OF WILSON CHAIN 

For definiteness, we consider the single-impurity An- 
derson model. It describes a spinful fermionic impurity 
level with energy and double occupancy cost U (with 
associated creation operators /q^, where /i =i, t denotes 
spin) , which acquires a level width T due to being coupled 
to a spinful fermionic bath with bandwidth W — I. Since 
the questions studied in this paper are of a generic na- 
ture and do not depend much on the specific parameter 
values used, we consider only the symmetric Anderson 
model and take U = ^, U/ttT = 1.013 and — —\U 
throughout this paper. To achieve a separation of energy 
scales, following Wilsonjii^ the bath is represented by a 
set of discrete energy levels with logarithmically spaced 
energies A~" (with associated creation operators /^^), 
where n>l,A>lisa "discretization parameter", and 
the limit A 1 reproduces a continuous bath spectrum. 
The discretized Anderson model Hamiltonian can then 
be represented as 

Ham = lim Hn , (1) 

N^oo 

where Ti.^ describes a Wilson chain of "length N" (i.e., 
up to and including site N): 

Hn = T^Ni+Ti-Ni +U{f^^ff^^f^^ff^^ + (2a) 

with hopping coefficients given by 

t^^i^ f°'^" = 0' (3) 

ii(l + A-i)A-("-i)/2e, forn>l, 

^„ = (1- A-")(l~A-2»+i)-i/2(^_^-2„-i)-i/2 ^ 

In passing, we note that for our numerics we have found 
it convenient (following Refs. [l3 and [iB) to keep track 
of fermionic minus signs by making a Jordan- Wigner 
transformationist of the Wilson chain to a spin chain, us- 
i^S fnfj. ~ ^n^iS^fj, and fnf_i = Pn^iSj^^. Here are a set 
of spin-i raising and lowering operators, that for equal 
indices satisfy {s-^,s+^} = 1, {s'^)^ = [s+^f = 0, but 
commute if their indices are unequal. The fermionic an- 
ticommutation relations for the fn^ are ensured by the 
operators = (— l)^(fiA)<(nM) "ftp^ftA^ where < refers to 
some implicitly specified ordering for the composite index 
{njj) . The need to be kept track of when calculating 
certain correlation functions, but do not arise explicitly 
in the construction of the matrix product states that are 
the focus of this paper. This transformation will implic- 
itly be assumed to have been implemented throughout 
the ensuing discussion. 

For the Anderson model, site n of the Wilson chain rep- 
resents the set of four states |(t„), with cr„ = (cnXi'^nT) ^ 



{(00), (10), (01), (11)}, where cr„,, e {0, 1}, to be viewed 
as eigenvalue of s+^s^^, gives the occupancy on site n 
of electrons with spin /x. Thus, the dimension of the 
spinful index cr„ is = 4, and that of the spin-resolved 
index cr„^ is d! = 2. As a general rule, we shall use the 
absence or presence of primes, d vs. d' (and D vs. D' 
below), to distinguish dimensions referring to spinful or 
spin-resolved indices, respectively, and correspondingly 
to folded or unfolded representations of the Wilson chain. 
For other quantum impurity models, such as the Kondo 
model or multilevel Anderson models, the dimension of 
the local impurity site, say do, differs from that of the 
bath sites, do ^ d. It is straightforward to generalize the 
discussion below accordingly. 

The Hamiltonian Tijv of a Wilson chain of length N 
is defined on a Hilbert space of dimension d^^^. It is 
spanned by an orthonormal set of states that, writing 
|cn) = \^ni)\<^n^) , Can be written in either spinful or 
spin-resolved form, 

|^^> - |ao)|ai)...|a^), (4a) 

= koi)|CToT)|CTu)|o-iT)...|o-Ari)|crAr|), (4b) 

corresponding to a "folded" or "unfolded" representation 
of the Wilson chain, illustrated by Figs. [DJa) or (b), 
respectively. The unfolded representation of Fig. [T]Jb) 
makes explicit that the Anderson Hamiltonian of Eq. ([2]) 
has the form of two separate Wilson chains of specified 
spin, described by TIni ^'^d H-n^, which interact only 
at site zero. This fact will be exploited extensively be- 
low. Note that the ordering chosen for the |(T„^) states 
in Eq. ijib)) fixes the structure of the many-body Hilbert 
space once and for all. The fact that the sites of the 
unfolded chain in Fig. [T] are connected in a different or- 
der than that specified in Eq. (|4bp is a statement about 
the dynamics of the model and of no consequence at this 
stage, where we simply fix a basis. 



III. NRG TREATMENT OF FOLDED WILSON 
CHAIN 

A. NRG matrix product state arises by iteration 

Wilson proposed to diagonalize the folded Wilson chain 
numerically using an iterative procedure, starting from 
a short chain and adding one site at a time. Consider 
a chain of length rt, sufficiently short that Tin can be 
diagonalized exactly numerically. Denote its eigenstates 
by |i?2)f, ordered by increasing energy {E2)[, with a — 
1, . . . , D„ and Z?„ = c?"+^. (We use subscripts f and u to 
distinguish quantities obtained from a folded or unfolded 
Wilson chain, respectively; similarly, in later parts of the 
paper we will use the subscripts r and c for "refolded" 
and "cloned".) E.g., for a chain consisting of only the 
impurity site, n = 0, the d eigenstates can be written as 

linear combinations of the form \E^){ = J^ao \'^o)A^ia\ 
where the coefficients have been arranged into d matrices 
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(a) 




Figure 1: (Color online) (a) The standard spinful or "folded" 
representation of the Wilson chain of the single-impurity An- 
derson model, and (b) its spin-resolved or "unfolded" repre- 
sentation. The latter makes explicit that spin-down and -up 
states are coupled only at the impurity sites and not at any 
of the bath sites. The dashed boxes indicate the chains de- 
scribed by Hi and Hn, respectively. 

^["■o] Qf dimensions 1 x d (i.e., d-dimensional vectors), 
with matrix elements A^[°\ Then add to the chain the 
site n + 1 and diagonalize Hn+i in the enlarged Hilbert 
space spanned by the (Dnd) states |i?2)t|cr„+i). The new 
orthonormal set of eigenstates, with energies {Ep~^^){, 
can be written as linear combinations of the form 

d Dr, 
— 1 a—1 

with l3 = 1, . . . , (Dnd). Here the coefficients specifying 
the linear combination have been arranged into a set of 
d matrices of dimension Z3„ x Dn+i, with matrix 

elements A^^^^^^^- The orthonormality of the eigenstates 
at each stage of the iteration, f{EJ^\Ep,)f = Sppi, implies 
that the A-matrices automatically satisfy the orthonor- 
mality condition 

^^K]t^k„]^j_ (6) 

We remark that it is possible to exploit symmetries of 
Tin (e.g. under particle-hole transformation) to cast A in 
block-diagonal form to make the calculation more time- 
and memory-efficient. However, for the purposes of the 
present paper, this was not required. 

Iterating the above procedure by adding site after site 
and repeatedly using Eq. |[5|), we readily find that the 



NRG eigenstates of T-Ln on the folded Wilson chain can be 
written in the form of a so-called matrix product stated 

\E^)i = E l*^"^) (^'"'''^'"^' ■ • • , (7) 

illustrated in Fig. [2ta). Here matrix multiplica- 
tion is impHed in the product, {A^'^"-^ A^""^+^^)ap — 

Aaj^ aI!^^'''^\ and {cr^} denotes the set of all se- 
quences ctq, (Ti, . . . , (Tat. This matrix multiplication gen- 
erates entanglement between neighboring sites, with the 
capacity for entanglement increasing with the dimension 
Dn of the index being summed over. 



B. NRG truncation 

In practice, it is of course not possible to carry out the 
above iteration strategy explicitly for chains longer than 
a few sites, because the size of the ^-matrices grows ex- 
ponentially with N. Hence Wilson proposed the follow- 
ing NRG truncation procedure: Once Dn becomes larger 
than a specified value, say D, only the lowest D eigen- 
states \E2){, with a — 1, . . . , D, are retained or kept at 
each iteration, and all higher-lying ones discardedi^. Ex- 
plicitly, the upper Hmit for the sum over a in Eq. ^ is 
redefined to be 

£>„ = min(d"+\ £») . (8) 

As a result, the dimensions of the A^"^"' matrices occur- 
ring in the matrix product state ^ start from 1 x d at 
n = and grow by a factor of d for each new site until 
they saturate at D x D after truncation has set in. The 
structure of the resulting states \E^)f is schematically 
depicted in Figs.[2l[a) and[2l[b), in which the site index is 
viewed as a single or composite index, (t„ or (cr„|,cr„|), 
respectively. 

Wilson showed that this truncation procedure works 
well in practice, because the hopping parameters tn of 
Eq. (IH) decrease exponentially with n: the resulting sep- 
aration of energy scales along the chain ensures that high- 
lying eigenstates from iteration n make only a small con- 
tribution to the low-lying eigenstates of iteration n+1, so 
that discarding the former hardly affects the latter. The 
output of the NRG algorithm is a set of eigenstates \EJ^)f 
and eigenenergies (i?^)f for each iteration, describing the 
physics at energy scale A~"/^. The NRG eigenenergies 
are usually plotted in rescaled form, 

(e^)f = {E^ - i?r)f/A-"/2 , (9) 

as functions of n, to obtain a so-called NRG fiow diagram; 
it converges to a set of fixed-point values as n oo. 
Figure[7]in Section lV BI below shows some examples. The 
ground state energy of the entire chain is given by the 
lowest energy of the last iteration, {Eq)i = {E^){. 

Despite the great success of NRG, Wilsonian trunca- 
tion is does have some drawbacks. Firstly, its errors grow 
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Figure 2: (Color online) (a) and (b) show the matrix prod- 
uct structure of the state \E^)f of Eq. ([T]), depicting the site 
index as a single or composite index, a„ or (cr„|, (J„|), respec- 
tively, (c) shows the matrix product structure of the state 
|*^)u of Eq. dl). (For the sake of illustrating Eq. of 
Appendix IA 2\ the labels {Bni)vi,i in the bottom row are pur- 
posefully typeset "upside down", so that they would be right- 
side up if the chain of boxes were all drawn in one row in the 
order indicated by Eq. ifTSjl .l Each matrix yl or B is repre- 
sented by a box, summed-over indices by links, free indices 
by terminals, and dummy indices having just a single value, 
namely 1, by ending in a triangle. The dimensions (d, D, 
d , D , etc.) next to each link or terminal give the number of 
possible values taken on by the corresponding index, assuming 
Wilsonian truncation for (a) and (b), and VMPS truncation 
for (c). Note the similarity in structure between (c) and (b): 
the dashed boxes in the former, containing ® ^!,''^^' i 

play the role of the Aj^'f^ ^ matrices in the latter. Their 
capacity for entangling neighboring sites is comparable if one 
chooses -D'^ oc D [cf. Eq. ifM]! ]. since neighboring dashed boxes 
in (c) are connected by two links of combined dimension D'^, 
whereas neighboring A-matrices in (b) are connected by only 
a single link of dimension D. 



systematically as A tends to 1, because then the separa- 
tion of energy scales on which it relies becomes less effi- 
cient. Secondly, it is not variational, and hence it is not 
guaranteed to produce the best possible approximation 
for the ground state within the space of all matrix prod- 
uct states of similar form and size. We shall return to 
this point later in Section [Vl] and study quantitatively to 
what extent the NRG ground state wavefunction can be 



improved upon by further variational optimization. 



C. Mutual information of opposite spins on site n 

A crucial feature of the folded Wilson chain is that 
all degrees of freedom associated with the same energy 
scale, A~"/^, are represented by one and the same site 
and hence are all added during the same iteration step. 
Since the spin-down and -up degrees of freedom associ- 
ated with each site are thus treated on an equal footing, 
the resulting matrix product state provides comparable 
amounts of resources for encoding entanglement between 
local states of the same spin, involving |(T„^)|(T„+i^), or 
between states of opposite spin (indicated by the bar), 
involving |cr„^)|cr„p) or |cr„p) |cr„+ip) . However, it turns 
out that for the Anderson model this feature, though a 
priori attractive, is in fact an unnecessary (and memory- 
costly) luxury: Since the Anderson model Hamiltonian 
(121) couples spin-down and -up electrons only at the im- 
purity site, the amount of entanglement between states 
of opposite spin rapidly decreases with n. 

To illustrate and quantify this claim, it is instructive 
to calculate the so-called mutual information M^^ of the 
spin-down and -up degrees of freedom of a given site 
n. This quantity is defined via the following general 
construction Let C denote an arbitrary set of degrees 
of freedom of the Wilson chain, represented by the states 
\a^). Let be the reduced density matrix obtained 
from the ground state density matrix by tracing out all 
degrees of freedom except those of C, denoted by N\C: 



{<T«\C} 



N\Cx 



(10) 



For example, if C represents the spin-down and up- 
degrees of freedom of site n, its matrix elements are: 

pC^, ^ J2 (A['^"1^..^[^"1^..^[^°1^)gi 

x(A['"»1 ...v4f'"~l)iG . (11) 

If C represents only the spin-/i degree of freedom of site 
n, a similar expression holds, with n replaced by n/i. The 
entropy associated with such a density matrix is given by 



5^ 



Ini 



(12) 



= 1. 



where wf are the eigenvalues of p*", with 
Now, consider the case that C = AB is a combination of 
the degrees of freedom of two distinct subsets A and B, 
represented by states of the form = \a^)\a^). Then 
the mutual information of A and B, defined by 



M 



AB 



s 



AB 



(13) 



characterizes the information contained in p beyond 
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Figure 3: NRG result for the mutual information M^^ be- 
tween spin-down and -up degrees of freedom of site n of a 
folded Wilson chain of lenght N = 50. The Anderson model 
parameters are fixed at U = ^, U/nF — 1.013, €d — —\U 
throughout this paper. Lines connecting data points are 
guides for the eye. The slight differences in behavior observed 
for even or odd n are reminiscent of the well-known fact^ that 
the ground state degeneracy of a Wilson chain is different for 
even or odd A*'. 



j^fAB — Q jf there is no entanglement between the degrees 
of freedom of A and B, since then p^^ — ^ and 
its eigenvalues have a product structure, wf,^ = w{^w" 
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We define the mutual information between spin-down 
and -up degrees of freedom of site n of the folded chain, 
M^T, by Eq. taking A = ni and B = n]. Figure[l 
shows this quantity as function of n for the symmetric 
Anderson model. Evidently is very small for all but 
the first few sites, and decreases exponentially with n. 
This implies that for most of the folded chain, there is 
practically no entanglement between the spin-down and 
-up degrees of freedom. Consequently, the corresponding 
matrices occuring in Eq. ([7| for \EQ)i in effect have a 
direct product structure: loosely speaking, we may write 
^k„] ^ ^Ki] ^ In the next subsection, we will 

exploit this fact to achieve a significant reduction in mem- 
ory cost, by implementing the effective factorization in an 
alternative matrix product Ansatz [see Eq. ifTSi) below], 
defined on an unfolded Wilson chain which represents n | 
and nt of freedom by two separate sites. 



IV. DMRG TREATMENT OF UNFOLDED 
WILSON CHAIN 



A. Variational matrix product state Ansatz 

As pointed out by Verstraete et al.,^^ an alternative 
approach for finding a numerical approximation for the 
ground state can be obtained by variationally minimizing 
the ground state energy in the space of all "variational 
matrix product states" (VMPS) of fixed norm. Imple- 
menting the latter constraint via a Lagrange multiplier 
A, one thus considers the following minimization prob- 
lem. 



min [(^-IWatI*) 

l*>e{|*">u} 



(14) 



The minimization is to be performed over the space of all 
variational matrix product states |^'^)u having a speci- 
fied structure (see below), with specified dimensions 
for the matrices, whose matrix elements are now treated 
as variational parameters. This minimization can be per- 
formed by a "sweeping procedure", which optimizes one 
matrix at a time while keeping all others fixed, then opti- 
mizing the neighboring matrix, and so forth, until conver- 
gence is achieved. The resulting algorithm is equivalent 
to a single-site DMRG treatment of the Wilson chain. 
Our main goal is to analyse how the energies and eigen- 
states so obtained compare to those produced by NRG. 

Having decided to use a variational approach, it be- 
comes possible to explore matrix product states having 
different, possibly more memory-efficient structures than 
those of Eq. ^ and Fig. ^a). In particular, we can 
exploi1»i2, the fact that the Anderson model Hamiltonian 
([2]) couples spin-down and -up electrons only at the im- 
purity site, as emphasized in Eq. ^ and Fig. [l]Jb). For 
such a geometry, it is natural to consider matrix product 
states defined on the unfolded Wilson chain (subscript u) 
and having the following form, depicted schematically in 
Fig.EJc): 

= H |cr^)(S[''"il...5['"oils['"'"]...5['""Tl)^^. 

(15) 

The order in which the bI'^"^'! matrices occur in the prod- 
uct mimics the order in which the sites are connected in 
the unfolded Wilson chain. (The fact that this order 
differs from the order in which the basis states |(t„^) for 
each site are arranged in the many-body basis state |(t^), 
see Eq. ijib)) . does not cause minus signs compHcations, 
because we work with Jordan- Wigner-transformed effec- 
tive spin chains.) Each i?!'^"*'! stands for a set of d' = 2 

matrices with matrix elements i?!,'^"''', with dimensions 
Z?; X for Bl'^-il and D'^^_^ x D'^ for bI'^-tI, where 



Dl^ = min{d'''-'\D'), 



(16) 



as indicated on the links connecting the squares in 
Fig. [2l[c). This choice of matrix dimension allows the 
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outermost few sites at both ends of the unfolded chain to 
be described exactly (similarly to the first few sites of the 
folded Wilson chain for NRG) , while introducing trunca- 
tion, governed by D' , for the matrices in the central part 
of the chain. The first index on B^"J'^'^ and the second 
index on B^^^^^ are dummy indices taking on just a sin- 
gle value, namely 1, since they represent the ends of the 
chain. The triangles in Fig. [2l[c) are meant to represent 
this fact. As a result, Eq. (flSl) represents just a single 
state, namely the ground state, in contrast to Eq. 
which represents a set of states, labeled by the index /?. 
Moving inward from the endpoints by decreasing n, the 
matrix dimension parameter D'^ increases by one factor 
of d' for each site, in such a way that the resulting matri- 
ces are of just the right size to describe the outside ends 
of the chain (from n to TV) exactly^ i.e. without trunca- 
tion. After a few sites, however, truncation sets in and 
the matrix dimensions saturate at D' x D' for the central 
part of the chain. 

To initialize the variational search for optimal B- 
matrices, it turns out to be sufficient to start with a 
set of random matrices with normally distributed ran- 
dom matrix elements. Next, singular value decomposi- 
tion is used to orthonormaHze the B-matrices in such a 
way [see Eq. l|Aip ] that the matrix product state Eq. ifTSj) 
has norm 1 (see App. [XT] for details). Thereafter, vari- 
ational optimization sweeps are performed to minimize 
Eq. (fT4l) one i?-matrix at a timei^-. (The technical de- 
tails of this procedure will be published separately^.) 
After a sweeping back and forth through the entire chain 
a few times, the variational state typically converges (as 
illustrated by Fig. [13] in Sec. IVIB] below), provided that 
D' is sufficiently large. We shall denote the resulting con- 
verged variational ground state by \Eq)^. Its variational 
energy, {Eq)u, turns out to be essentially independent of 
the random choice of initial matrices. 



B. VMPS truncation 

Since D' x D' is the maximal dimension of iJ-matrices, 
D' is the truncation parameter determining the effective 
size of the variational space to be searched and hence the 
accuracy of the results. Its role can be understood more 
expHcitly using a technique that is exceedingly useful in 
the VMPS approach, namely singular value decomposi- 
tion: any rectangular matrix B of dimension m x m' can 
be written as 



B = USV'' , with U^U = Vt V = 1 



(17) 



where 5 is a diagonal matrix of dimension min(m, m'), 
whose diagonal elements, the so-called "singular values", 
can always be chosen to be real and non-negative, and 
U and are column- and row-unitary matrices (with 
dimensions m x min(m,m') and min(m, m') x to', re- 
spectively). Due to the latter fact, the matrix norm of B 
is governed by the magnitude of the singular values. 
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Figure 4: (Color online) (a) Typical singular value spectrum 
for site 5^ of the unfolded Wilson chain, obtained by singular 
value decomposition of B^°'^i\ It shows, roughly, power-law 
decrease for large enough /3, modulo steps due to degeneracies 
in the singular value spectrum, (b) D'-dependence of the 
truncation error t(D') [Eq. I|18p ]. 



For any given site of the unfolded Wilson chain, this de- 
composition can be applied in one of two ways (depending 
on the context, see App.[A]) to the set of matrices with el- 
ements siy^"^': introduce a composite index u — {anfi, v) 
(or fj = (ct„^,77)) to arrange their matrix elements into a 
rectangular matrix carrying only two labels, with matrix 

elements B^,^ = (or B^fi — -Si'r)"^')) cind decompose 

this new matrix as B — USVK 

Now, if this is done for any site for which the set of ma- 
trices i?!'^"''! have maximal dimensions D' xD', the corre- 
sponding matrix S will likewise have dimensions D' x D' . 
Let its diagonal elements, the singular values (with 
v = 1, . . . ,£)'), be labelled in order of decreasing size. 
(Their squares, si, correspond to the eigenvalues of the 
density matrix constructed in the course of the single-site 
DMRG algorithm^.) If D' is sufficiently large, the Si, are 
typically found to decrease with increasing v roughly as 
some negative power of i^, as illustrated in Fig.[4Ua). The 
last and smallest of the singular values, s'jj, (squared, fol- 
lowing Ref . , thus indicates the weight of the informa- 
tion that is lost at that site due to the given (finite) choice 
of D': by choosing D' larger, less information would be 
lost since more singular values (though of smaller size) 
would be retained. Repeating such an analysis for all 
sites of the unfolded Wilson chain, one may define the 
largest of the sjj, parameters of the entire chain. 



r{D') 



max(s^/) 

{rifj.} 



(18) 



as "truncation error" characterizing the maximal infor- 
mation loss for a given value of D' . Typically, t{D') 
decreases as some negative power of D' , as illustrated 
in Fig. \Mi>)- In this way, D' assumes the role of a cut- 
off parameter that directly governs the accuracy of the 
VMPS approach, in a way analogous to the parameter D 
of NRG. 
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C. Refolding 

The VMPS approach purposefully focusses on finding 
an optimal description of the variational ground state 
\Eq)u- Nevertheless, the B-matrices from which the lat- 
ter is constructed contain information about all energy 
scales of the model, due to the logarithmic discretization 
of the Wilson chain. In particular, information about the 
scale A~"/^ is encoded in the set of matrices i?!'^""! as- 
sociated with the two site rt| and n\. From these, it is 
possible to extract excited-state eigenspectra and energy 
flow diagrams in complete analogy to those produced by 
NRG. In this subsection we explain how this can be ac- 
complished by a technique to be called "refolding", which 
combines the two matrices i?!'^"^] and ijI'^^Tl into a single 
matrix, say B^''"\ and thereby recasts unfolded matrix 
product states into folded ones. It should be emphasized 
that this procedure simply amounts to an internal re- 
organization of the representation of the VMPS ground 
state. 

Consider a given matrix product state |^^)u of the 
form lfT5|) , deflned on an unfolded Wilson chain of length 
A'' (e.g. the converged ground state \Eq)u). To refold it 
(subscript r) , it is expressed as a state of the following 
form [same as Eq. ([7])] 

l*^>r= J2 , (19) 

deflned on a folded Wilson chain of length TV and nor- 
maHzed to unity, i.(^"|^")r = 1- Graphically speaking, 
this corresponds to rewriting a state of the form shown in 
Fig. (He) in terms of states of the form of Fig.jS^a). To ob- 
tain the matrices needed for Eq. lfT9|). one constructs, for 
every site n of the refolded chain, a set of d matrices Sl"'"! 
from a combination of the two sets of spin-resolved ma- 
trices ^['^"il and ^['^"t] of the unfolded chain (App. [Al 
gives the details of this construction). This is done in 
such a way, using singular value decomposition, that (i) 
the resulting matrices B^°'"^ satisfy the orthonormality 
conditions ^ (with A ^ B), thereby guaranteeing the 
unit normalization of the the refolded state l^*^)!; and 
(ii) the ijl"'"] matrices have a structure similar to that of 
the matrices A^'^^^ generated by NRG, except that their 
dimensions, D^^ x D^^^i, are governed by 

= min(d", d^+i-", D'^) (20) 

[instead of Eq. for reasons explained in App. lA 21 
Thus, their dimensions have the maximal value D'^ x D^, 
with = D'^, in the central part of the refolded chain, 
while decreasing at its ends towards 1 x d or c? x 1 for 
n = or iV, respectively. 

Now, suppose that a converged variational ground 
state \Eq)u has been obtained and refolded into the form 
|5'^)r, so that the corresponding orthonormaHzed matri- 
ces i?!'^"! for the refolded Wilson chain of length TV are the 
building blocks of the ground state of the system. Then 



it is possible to extract from them information about the 
many-body excitation spectrum at energy scale A~"/^ 
that is analogous to the information produced by NRG. 
To this end, consider a subchain of length n of the full 
refolded Wilson chain, and use the deflnition 

\^}),= J2 k")(i?["°ii?["^i...i?[""i)i/3 , (21) 

[as in Eq. (fT9|) . but with N replaced by n] to construct a 
set of states l^*^)!- on this subchain. These states, shown 
schematically by sites to n of Fig. [2l[a), form an or- 
thonormal set, r(^Q|^^)r = <5q^, due to the orthonor- 
mality [Eq. ^] of their constituent matrices. They can 
thus be viewed as a basis for that subspace of the many- 
body Hilbert space for the length-rt Wilson chain, i.e. of 
that subspace of span{|<T")}, which VMPS sweeping has 
singled out to be most relevant for describing the ground 
state \Eq)u of the full chain of length N. Therefore we 
shall henceforth call the |*;^)r "(refolded) VMPS basis 
states" for this subchain. 

This basis can be used to deflne an effective "refolded 
Hamiltonian" H" for this subchain, with matrix elements 

{Hi%p ^ An\nn\n)r ■ (22) 

Its eigenvalues and eigenstates, say (£'^')r and |i?^)r, are 
the VMPS analogues of the NRG eigenvalues and eigen- 
states, {EJ^)i and \EJ!j)[, respectively. They differ, in gen- 
eral, because VMPS and NRG use different truncation 
criteria, but are expected to agree well for sufflciently 
large choices of D' and D. This is indeed found to be the 
will be shown in detail in the next section. 



V. COMPARISON OF NRG AND VMPS 
RESULTS 

Having outlined the NRG and VMPS strategies in the 
previous section, we now turn to a comparison of their 
results. This will be done, in successive subsections, by 
comparing their ground state energies and the overlaps 
of the corresponding ground states; the eigenspectra and 
density of states obtained from both approaches; and fl- 
nally, the energy eigenstates used in the two approaches. 
We will thereby gain more insights into the differences 
between NRG and VMPS truncation criteria. Before em- 
barking on a detailed comparison, though, some remarks 
on the choices to be made for D and D' are in order. 

Since the structure of the matrix products occurring 
in Eqs. ^ and lfT5|) differ, the spaces consisting of all 
states of the type \EJ^)f or \E^)r, to be called the "NRG- 
subspace" or "VMPS-subspace" for a length-n chain, re- 
spectively, constitute nonidentical subspaces of the 
dimensional Hilbert space spanned by the basis states 
|<t"). The extent to which they describe the energy eigen- 
states oiH-N with comparable accuracy will depend very 
strongly on the choices made for D and D' . It turns out 
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(numerical evidence will be presented below) that with 
the choice 

D' = d'""^/D , (23) 

the VMPS-subspace is sufficiently large to give highly 
accurate representations of all kept states of NRG (in- 
cluding, in particular, the ground state) for the choice 
m = 1, or of all kept and discarded states of NRG for the 
choice TO = 2. The fact that D' should be proportional 
to ^/D can be made plausible by considering the follow- 
ing question: given a folded Wilson subchain of length n 
(i.e. consisting of sites to n) and its equivalent unfolded 
version, what are the smallest values for the dimensions 
D and D' for which both approaches describe the ground 
state exactly, i.e. without any truncation? Answer: On 
the one hand, the folded subchain has n + 1 sites of di- 
mension d, and hence a total dimension d""*"^; to ensure 
that the ground state in this space is described exactly, 
the kept space of the previous iteration must not involve 
any truncation, implying D = d^. On the other hand, for 
the equivalent unfolded subchain, the spin J, and t parts 
each have n + 1 sites of dimension d' , hence each have a 
Hilbert space of total dimension d'*-"^^"*; to ensure that 
this space is described without truncation, its dimension 
should equal the maximal dimension of the i?-matrices 
at sites 0/i, implying D' = d'"'*'^. Using d' = Vd we 
readily find D' — d'\/D, establishing the proportionality 
between D' and \/D and suggesting the choice m = 1 to 
achieve an accurate VMPS-representation of the ground 
state. Actually, we find numerically that already m = 
yields good qualitative agreement between the VMPS 
and NRG ground states, while to = 1 yields a quantita- 
tively accurate VMPS-representation of the NRG ground 
state also for larger chain lengths, that do involve trun- 
cation. Since such ground states are built from the kept 
spaces of previous iterations, this implies that for to = 1, 
all kept states in NRG (not only the ground state) are 
likewise well represented by VMPS. Indeed, we will find 
this to be the case. Moreover, it turns out numerically 
that with TO = 2, it is also possible to achieve an ac- 
curate VMPS-representation of all kept and discarded 
NRG-type states, as will be extensively illustrated be- 
low. 

For the results reported below, we show data only for 
even iteration number n, to avoid even/odd oscillation ef- 
fects that are typical and well-understood for Wilsonian 
logarithmic discretization, but not of particular interest 
here. We set D' = d'"^^fD throughout and specify the 
choices made for to. All VMPS results shown in this sec- 
tion are extracted from randomly initialized, fully con- 
verged variational ground states IS'g )u of the form (flSl) . 



A. Ground state energies and overlaps 

Figures EJa) and [Sljb) compare the NRG and VMPS 
ground state energies, {Eq)[ and (£'q )u, for three values 



of A and, in (a), two values of m. They illustrate three 
points. Firstly, for a given A the VMPS ground state 
energies are smaller than those of NRG, (-Eq )u < (-E-q )f j 
as expected, since VMPS is a variational method and 
NRG is not. Secondly, Fig. [H^a) shows that larger val- 
ues of TO yield lower VMPS ground state energies, as 
expected, since their variational space is larger. Thirdly, 
the improvement of VMPS over NRG, as measured by 
the energy difference (-Eq )f — {Eq)u shown in Fig.[5l[b), 
becomes more significant for smaller A, as expected, since 
the truncation scheme of NRG relies heavily on energy 
scale separation, and hence becomes less efficient for 
smaller A. 

Figure [DJc) compares the overlap between NRG and 
VMPS ground states, characterized by the deviation from 
1 of the overlap |t(i?Q )u|. The latter can be calcu- 
lated straightforwardly from 

t{E^ \E^),, = J2 (^'""'^ ■ • • ^'""'^)gi (24) 
X (i?['^"il . . . ijl'^oil ^kml ^ ^ ^ sI'^wtI)^^ 

where the index contractions associated with the sum- 
mation over repeated indices are illustrated in Fig. [6Ua). 
Fig. (He) shows that the deviation of the overlap from 1 
becomes larger the smaller A, again illustrating that then 
the NRG truncation scheme becomes less reliable. 



B. Comparison of eigenspectra and density of 
states 

Figure [7] compares the energy flow diagrams obtained 
from NRG and refolded VMPS data, the latter obtained 
by diagonalizing the effective Hamiltonian of Eq. (|22l) . 
It shows the rescaled energies (e^)f,r of Eq. ([9]) as func- 
tions of n, for four combinations of m and A, and illus- 
trates the same trends as found in the previous subsec- 
tion: Firstly, the NRG and VMPS fiow diagrams clearly 
agree not only for the ground state but also for a signifi- 
cant number of excited states. Evidently, the variational 
space searched by VMPS is large enough to capture con- 
siderable information about excited states, too, although 
the VMPS method was designed to optimize only the 
ground state. Moreover, for a given choice of A, NRG 
and VMPS eigenenergies coincide for a larger number of 
states for m = 2 than for to = [compare (b) to (a) and 
(d) to (c)], because the variational space is larger. Sec- 
ondly, for a given choice of m, NRG and VMPS eigenen- 
ergies agree better for A = 2.5 than for A = 1.5 [compare 
(c) to (a) and (d) to (b)], as expected, because larger 
A leads to better energy scale separation and reduces 
the inaccuracies inherent in NRG's Wilsonian truncation 
scheme. 

As a complementary way of analysing spectral infor- 
mation we also consider the "density of states", for a given 
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Figure 5: (Color online) Comparison of NRG and VMPS re- 
sults for (a,b) the ground state energies and (c) the ground 
state overlaps, plotted as a functions of D with D' = d'"^y/D, 
for three values of A and, in (a), for two values of m. In (a) the 
reference energies for each A were obtained by extrapo- 
lating the VMPS data points for m = 2 to D' ^ cxa, which 
represents the best estimate of the true ground state energy 
available within the present set of methods. The power law 
fits to the numerical data in (b) and (c), shown as dashed 
lines, were made for the three data points with largest D, for 
which the dimensions are large enough to have reliable NRG 
data. 
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Figure 6: (Color online) Contraction patterns used to cal- 
culate (a) the overlap f{-BG I-E'g)u [Eq. lf24|l ] between folded 
NRG and unfolded VMPS ground states, and (b) the overlap 
matrix = r(*S|-E^')t [Eq. ifm]! ] between refolded VMPS 
basis states and folded NRG eigenstates. Boxes represent 
A or B matrices in the graphical representation of Fig. [2] 
and links connecting them represent indices that are being 
summed over. 



-la- 



I [a\pK) is a Gaussian peak of width a 



using the rescaled eigenenergies of Eq. ([9]). Here 



and unit weight, used to broaden the discrete spectrum 
in order to be able to plot it, and the number of states 
included in the sum is taken as Umax — dD or d"^D for 
NRG or VMPS results, respectively. Figure [8] shows such 
a density of states for several choices of m and iteration 
number n. It illustrates three points: 

Firstly, although for small energies p„(e) grows rapidly 
with e, as expected for a many-body density of states, it 
does not continue to do so for larger e (the exact density 
of states would) , due to the truncation inherent in both 
NRG and VMPS strategies. For NRG, p„(e) drops to 
very abruptly, because by construction Wilsonian trun- 
cation is sharp in energy space (at each iteration only 
the lowest dD eigenstates are calculated). In contrast, 
for VMPS Pn(e) decreases more gradually for large e, 
because VMPS truncation for states at site n is based 
not on their energy, but on the variationally determined 
weight of their contribution to the ground state of the 
full Wilson chain of length N. Evidently, these weights 
decrease with increasing e less rapidly than assumed by 
NRG. 

Secondly, the agreement of the VMPS curve for Pni^) 
with that of NRG is rather poor for m = (disagreement 
sets in already within the range of kept states of NRG, 
indicated by the shaded region), better for m = 1 (the 
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Figure 7: (Color online) Comparison of energy flow diagrams from NRG (dashed red lines) and refolded VMPS data (solid 
black lines), showing the reseated energies (e^)f,r [Eq. ((O}] versus n, calculated for even iteration numbers and four combinations 
of m (= or 2) and A (= 1.5 or 2.5). The number of NRG states shown (kept and discarded) is Dd\ the number of refolded 
VMPS states shown is D' = D'^ = dTD, this being the maximal dimension of refolded matrices S'"^"'. For m = 2 and A = 2.5, 
the NRG and DMRG flow diagrams agree very well, see (d). 



range of kept states is fully reproduced) , and very good 
for TO = 2 (disagreement sets in only close to the upper 
end of range of discarded states) . 

Thirdly, for large n, Pn{£) becomes increasingly spiky. 
This reflects the fact that the spectrum approaches a 
fixed point with regularly-spaced eigenenergies, as is ev- 
ident in the energy fiow diagrams of Fig. [7l 



C. Comparison of energy eigenstates 

To compare the energy eigenstates produced by NRG 
and refolded VMPS for a chain of length n, we analyse 
the overlap matrix 

Sl,^,{El\E-),. (26) 



It can be conveniently calculated from 5" = U^S^, 
where U^p = r{E'^\'^0)i- is the matrix that diagonalizes 
the effective Hamiltonian matrix H^p of Eq. (|22l) . and 
the matrix 

S-^ ^ ,{n\E-ph , (27a) 
= J2 (S['^"lt,..s['^°l^)ai(A["°l...A[""l)i;3 (27b) 

characterizes how much weight the NRG eigenstates have 
in the space spanned by the refolded VMPS basis states, 
and vice versa. The contractions implicit in Eq. (|27bp 
are illustrated in Fig. [HUb). 

Figure [9] shows the overlap matrix S^p on a color scale 
ranging from to 1, for m = 1 and several values of n. 
For the region of low excitation energies (about the first 
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Figure 8: (Color online) Results for the density of states, Pn(e) 
[Eq. l[25ll ]. broadened with a Gaussian broadening function. 
In each panel, the red vertical dashed and solid lines [which 
coincide in (a)] indicate the energies of the highest-lying kept 
and discarded NRG states of that iteration, while the shaded 
area indicates the range of kept NRG states. 



hundred or so states) its structure is evidently close to 
block-diagonal, indicating that both sets of states from 
which it is built are reasonably good energy eigenstates. 
Had both sets been perfect energy eigenstates, as would 
be the case for D' and D large enough to avoid all trun- 
cation, the blocks would be completely sharp, with sizes 
determined by the degeneracies of the corresponding en- 
ergies. Sharp blocks are indeed observed for n = 2 
[Fig- in a)], because no truncation has occurred yet. The 
"fuzziness" shown by the blocks in Fig. ^h) to [9l[d) for 
larger n implies that truncation is beginning to make it- 
self felt, causing NRG and VMPS to increasingly disagree 
on how to construct the eigenstates corresponding to a 
given range of eigenenergies. Note that the fuzziness be- 
comes markedly more pronounced for a, (3 > 256. The 
reason is that whenever S^p is nonzero for (3 > D, the as- 
sociated VMPS states have weight among the discarded 
states of NRG, implying that NRG discards some states 
relevant for building the VMPS ground state. Thus, S^p 
quite literally measures to what extent the truncation cri- 
teria of NRG and VMPS are compatible. Near the end of 
the chain, for n = 18 [Fig. [9l[d)], the off-diagonal spread 
is significantly reduced compared to the middle of the 
chain {n = 6,12) [Fig. [9l[b,c)], for two reasons. Firstly, 
the dimensions of the refolded S-matrices become small 
for n near TV, see Eq. (|20|) . so that the amount of trunca- 
tion is much less severe near the end of the chain than in 
its middle. Secondly, the eigenspectra have converged to 
their fixed point values, so that the number of different 
eigenenergies in a given energy interval is reduced, thus 
reducing the fuzziness in Fig.[9l[d). 

Next consider the total weight which a given NRG- 
state \EJ^){ has within the refolded VMPS-subspace for 
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Figure 9: (Color online) Plot of the overlap matrix = 
r{KI-E;3)f [Eq. (dSl)] between refolded VMPS and NRG en- 
ergy eigenstates, with a color scale ranging between and 1. 
In (a), with n = 2, no truncation occurs at all, and both state 
labels a and (3 run from 1 to d"^^ = 64. In (b) to (d), trun- 
cation does occur: For the folded NRG eigenstates |i?^)f, the 
label a runs from 1 to Dd = 1024, i.e. it includes all kept and 
discarded NRG states, while for the refolded VMPS eigen- 
states \EJ^)r, the label /3 runs from 1 to = D'^ = 1024 
[Eq. ^]. 
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(28) 



It satisfies < < 1. Weights less than 1 im- 

ply that the VMPS-subspace is too small to adequately 
represent the corresponding NRG state. The second 
equality in Eq. l(28|) . which follows from the unitarity 
of U, is useful since it implies that these weights can 
also be calculated directly from the refolded states l^^r 
[Eq. ([2T|) ]. without the need for diagonalizing the large 
{D'^ X Z3'^-dimensional) effective refolded Hamiltonian 
[Eq. mi 

Figure [TOl shows such weights wjj"^ for various choices 
of rt, A and m. Their dependence on m reinforces the con- 
clusions of the previous subsection: For to = (blue -I- 
symbols) , the weights are equal to 1 for the lowest state of 
each iteration, but less than 1 for many of the kept states. 
This shows that the VMPS subspace is large enough to 
accurately represents the NRG ground state, but signif- 
icantly too small to accurately represent all kept states. 
For m = 1 (green x symbols), the weights are close to 
1 only for the kept states, while smoothly decreasing to- 
wards for higher-lying discarded states. Finally, for 
TO = 2 (orange o symbols) , the weights of both kept and 
discarded NRG states are all close to 1, implying that the 
VMPS subspace is large enough to accurately represent 
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Figure 10: (Color online) For several NRG iteration numbers 
n and two values of A (different panels), this figure shows the 
weights w^"' [Eq. ^] with which NRG states \Ep)f with 
rescaled NRG eigenenergies (ejj)f [Eq. JO])] are found to lie in 
the VMPS-subspace of dimension D' = with m = 0, 

1 or 2 (indicated by +, x or o, respectively). In each panel, 
the red vertical dashed and solid lines indicate the energies 
of the highest-lying kept and discarded NRG states of that 
iteration. For n = 3, both of these lines are missing, since 
truncation has not yet set in. The choices for n in the left 
and right panels of each row are related by A^''^^^ = A^""/^ 
to ensure that both panels show a comparable energy scale. 



essentially all states kept track of by NRG. Note that for 
m = and 1, the decrease of the weights w^"^ with in- 
creasing energy occurs in a smooth and gradual fashion, 
illustrating yet again the smooth nature of VMPS trun- 
cation when viewed in energy space. When a smaller 
value of A is used [compare panels (a-d) to (e-h)] the 
weights of the higher- lying states of a given iteration tend 
to spread out over a larger range of values, since NRG 
has a weaker energy scale separation for smaller A. Fi- 
nally, the increasing spikyness of the eigenspectrum with 
increasing n, see Fig.fTOlfd.h). is due to the approach to a 
fixed point spectrum with regularly-spaced eigenenergies, 
as mentioned above. 




'0 8 16 24 32 40 48 4 
n (NRG iteration) 



8 12 16 20 24 28 
n (NRG iteration) 



Figure 11: (Color online) Integrated weights (see 
Eq. lf29|l l for two different A and three values of m. Dashed 
lines depict the maximum possible values of the kept and dis- 
carded weights, J and | for and WS, respectively. 



The results just discussed may be represented more 
compactly by considering, for a given iteration n, the 
integrated weights obtained by summing up the weights 
of all NRG states of type X, 



in) 



1 



/3ex 



w 



in) 



(29) 



where X — K,D,A stands for kept, discarded or all, re- 
spectively. All three types of integrated weights are nor- 
malized to the total number dD of all NRG states cal- 
culated at a given iteration (with d — A here), and reach 
their maximal values (-j, | and 1, respectively) when all 
the individual weights for that iteration equal 1. Fig- 
ure [TT] shows such integrated weights for several values 
of m and A. Upon increasing m from to 2, the inte- 
grated weights tend toward their maximal values, doing 
so more rapidly for larger A. For m = 2, they essentially 
saturate their maximal values, indicating yet again that 
the VMPS variational space is now large enough to fully 
retain all information kept track of by NRG. 

To summarize the result of this section: The VMPS 
approach reproduces NRG ground state properties much 
more cheaply, requiring only D' = ^/D for qualitative 
agreement, and D' = d'^/T) for quantitative agreement. 
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Moreover, it can also reproduce all kept and discarded 
NRG eigenstates if D' = d!'^\fD is used. However, 
to obtain excited energy eigenstates, we have to refold, 
requiring the diagonaHzation of matrices of dimension 
D'^ X D'^ . The numerical cost of doing so is comparable 
to that of NRG. 

The fact that VMPS gives access to the same informa- 
tion on eigenstates and eigenvalues as NRG has a very 
significant and reassuring consequence: all physical prop- 
erties of the model that can be calculated by NRG can 
also be calculated by VMPS, in combination with refold- 
ing. 



VI. CLONING AND VARIATIONAL 
IMPROVEMENT OF NRG GROUND STATE 




D' (number of states kept in VMPS) 



Viewed in MPS language, the NRG method constructs 
the ground state in a single sweep along the chain: each 
A is calculated only once, without allowing for possible 
feedback of information from ^'s describing lower ener- 
gies to those of higher energies calculated earlier. Thus, 
the resulting NRG ground state |i?G)f, to be denoted 
simply by |G)f below, is not optimal in a variational 
sense. In this section we investigate to what extent the 
ground state energy can be lowered further by perform- 
ing variational energy optimization sweeps on |G)f that 
serve to account for feedback of information from low to 
high energy scales. This feedback turns out to be small 
in practice, as will be seen below, but it is not strictly 
zero and its importance increases as the logarithmic dis- 
cretization is refined by taking A 1 . 



A. Mapping folded to unfolded states by cloning 

Our first step is to rewrite a given NRG ground state 
|G)f in a form amenable to subsequent energy minimiza- 
tion sweeps. To this end, we use a variational cl 
procedure (subscript c). 



|G)f='^'^|G)ee{|*^M 



(30) 



which maps |G)f of the form of Eq. ([7]) [Fig. EI a)] onto 
an unfolded state |G)c of the form |*^)u of Eq. ^ 
[Fig. [2ljc)]. Since their matrix product structures differ, 
this mapping will, for general values of D and Z?', not be 
exact, though its accuracy should improve systematically 
with increasing D' and hence increasing dimensions of the 
variational space. To be explicit, we seek the best pos- 
sible approximation to |G)f in the space of all unfolded 
states of the form ifTSj) , by solving the minimization prob- 
lem 

min [II |G)f-|G)c||'+A(|||G)c||'-l)] , (31) 

which minimizes the "distance" between |G)c and |G)t 
under the constraint, implemented using a Lagrange mul- 



Figure 12: (Color online) The deviation of the overlap 
|c{G|G)f I from 1 (red circles) and the cloning truncation error 
t(D') (blue squares), as functions of the number D' of kept 
states in the cloning procedure. Both approach in power-law 
fashion, as indicated by the dashed line fits. The inset shows 
how the overlap deviation from 1 decreases and converges to 
a small but finite constant in the course of sequential cloning 
sweeps. 



tiplier A, that the norm c(G|G)c — 1 remains constant. 
Varying Eq. I|3ip with respect to the matrices defining 
|G)c leads to a set of equations, one for each kfi, of the 
form 



d 



(l + A)c(G|G)c-2Re(f(G|G)c) =0, (32) 



which determine the i?-matrices of the desired "cloned" 
state |G)c. These equations can be solved in a fashion en- 
tirely analogous to energy optimization: Pick a particular 
site of the unfolded chain, say kfi, and solve the corre- 
sponding Eq. (|32|) for the matrix B^'^'"^^ while regarding 
the matrices of all other sites as fixed. Then move on to 
the neighboring site and in this fashion sweep back and 
forth through the chain until convergence is achieved. 
Appendix I A 31 describes some details of this procedure. 

A figure of merit for the success of cloning is the de- 
viation of the overlap |c(G|G)f| from 1. This deviation 
decreases monotonically with successive cloning sweeps 
and converges to a small but finite (C-dependent) value 
when the cloning process converges, as illustrated in the 
inset of Fig. [121 The main part of Fig. [12] shows that 
when the number D' of VMPS states is increased, the 
converged value of the overlap deviation approaches as 
a power law in D' (red circles) . It also shows that the cor- 
responding VMPS truncation error t(D') incurred during 
cloning (blue squares), calculated according to Eq. (fT8| . 
likewise decreases in power-law fashion with D' . All in 
all. Fig. [12] confirms that cloning works very well if D' is 
sufficiently large. 
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Figure 13: (Color online) Comparison of (a) energies and (b) 
wave function overlaps for random initialization (squares) vs. 
NRG-cloned initialization (circles), as functions of the num- 
ber fccm of variational energy minimization sweeps. Results 
are shown for A = 1.5 (green, open symbols, dashed lines) and 
A — 2.5, (blue, filled symbols, solid lines). The energies in (a) 
and overlaps in (b) are calculated with respect to a reference 
ground state |G)ref with D' = 64, obtained by performing 50 
energy minimization sweeps starting from random initializa- 
tion. The red horizontal straight lines in (a) (dashed or solid 
for A = 1.5 or 2.5, respectively), show the energy difference 
i^NRG — fi'rof, where Snrg is the energy of the NRG ground 
state |G)f used as input into cloning. The fact that i?NRG 
does not completely coincide with the energy Ec = -Efcom=o of 
the cloned state (horizontal straight lines do not meet circles 
at fcem = 0) is due to the fact that the deviation of the overlap 
|c(G|G)f| from 1 is not strictly equal to (see Fig. I12p . 



B. Variational energy minimization of |G)c 

Having used cloning to find the optimal unfolded rep- 
resentation |G)c of the NRG ground state |G)f, we 
now variationally minimize its energy by sweeping. We 
thereby obtain a sequence of states jG)^"" of ever lower 
energy, Ek^^, where the index kcm = 0, 1, 2, . . . gives the 
number of energy minimization sweeps that have been 
performed. The procedure is precisely analogous to that 
described in Section HV A\ the only difference being that 
the random initial state used there is here replaced by 
the cloned state |G)J? = |G)c as initial state. 

Figure [TSlfa) shows the evolution of the ground state 



energy E^^^^ as function of the number fccm of energy min- 
imization sweeps, for both random (squares) and cloned 
(circles) initial states. i?fc„„ is displayed with respect 
to the energy i^rcf of a reference state |G)rcf , defined in 
the figure caption, which represents our best approxima- 
tion to the true ground state. Figure flSlfb) shows how 
1 — |rcf(G|G)c| decreases as sweeping procedes, converg- 
ing to a small but finite value. For a given value of A (1.5, 
shown in green, open symbols connected by dashed lines, 
or 2.5, shown in blue, filled symbols connected by solid 
lines), the energies for random and cloned initialization 
shown in Fig. fTSF a) converge to the same value within 
just a few sweeps. However, the convergence is quicker 
for the cloned (circles) than the random (squares) in- 
put state, since the former represents an already rather 
good initial approximation (namely that of NRG) to the 
true ground state, whereas the latter is simply a random 
state. Nevertheless, the circles show strikingly that the 
NRG ground energy is not optimal, in that the energy 
can be lowered still further by sweeping. Moreover, this 
improvement is more significant for small than large A 
(for circled data points, compare dashed green to soHd 
blue lines for A = 1.5 or 2.5, respectively). The reason is 
that the NRG truncation scheme becomes less accurate 
the smaller A is, implying that the NRG result can be 
improved more significantly by further sweeping. This is 
again a reminder that the systematic error of NRG in- 
creases as A approaches 1, as already observed in Fig.[5l 



VII. CONCLUSIONS 

In this paper we presented a systematic comparison 
between NRG and DMRG, which we mainly referred to 
as VMPS, for the single- impurity Anderson model within 
the framework of matrix product states. We first refor- 
mulated both NRG and DMRG in the language of MPS, 
using a folded Wilson chain for NRG and an unfolded one 
for DMRG. Then we quantitatively compared the results 
of NRG and the VMPS approach for energy eigenvalues 
and eigenstates and explicitly analysed the difference in 
their truncation criteria, which are sharp or smooth in 
energy space, respectively. 

The most important conclusion of our study is this: For 
the purpose of obtaining the ground state of this model, 
the VMPS approach applied to the unfolded Wilson chain 
yields a very significant increase in numerical efficiency 
compared to NRG {D' — d'^/D), essentially without loss 
of relevant information. The physical reason is that the 
spin-down and -up chains are only weakly entangled for 
this model, so that the NRG matrices yll'^"! of dimension 
D that describe site n of the Wilson chain, can, in effect, 
be factorized as a direct product 3^"^"-^^ ® i?['^"T] of two 
matrices, each having dimension d'^/lj. It should be em- 
phasized, though, that this property relies on the physics 
of the model, namely the weak entanglement of the spin 
down and up chains. To what extent this property sur- 
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vives for other impurity models should be a subject for 
further research, the two-channel Kendo model being a 
particularly interesting candidate in this respect. 

Nevertheless, the possibility of using unfolded Wilson 
chains to reduce numerical costs for ground state cal- 
culations is very attractive for possible applications of 
the VMPS method to more complicated models involv- 
ing more than one conduction band.-^i For example, the 
conductance through a quantum dot coupled to two leads 
can under certain conditions (linear response, zero tem- 
perature, Fermi liquid behavior, etc.) be expressed in 
terms a set of phase shifts that are uniquely determined 
by the ground state occupation of the dot energy levels f?- 
Thus, in such situations reliable knowledge of the ground 
state is sufficient to calculate transport properties. 

Going beyond ground state properties, we showed that 
the entire excited state eigenspectrum of both kept and 
discarded NRG states can be recovered within the VMPS 
approach with at least the same accuracy as NRG, by us- 
ing D' = d'^ X ^/D and refolding. However, the latter 
step requires a subsequent additional diagonalization of 
matrices of dimensions D'^ , giving rise to a significant in- 
crease in numerical resources compared to the case that 
only ground state information is required. A quantita- 
tive comparison between NRG and VMPS for the eigen- 
spectrum's energies and eigenstates showed better agree- 
ment for A = 2.5 than 1.5, due to the fact that the NRG 
truncation scheme becomes increasingly less accurate the 
closer A approaches 1. 

Finally, we used a cloning procedure to recast a given 
folded NRG ground state into an unfolded form, and 
showed that its energy could be lowered further by sub- 
sequent energy minimization sweeps. As expected, we 
found that sweeping improves the relative accuracy with 
which the ground state energy can be determined, the 
more so the smaller the value of A. For example, for 
A — 1.5 the accuracy changed from O(10~'*) before 
sweeping to 0(10"'') thereafter [see Fig.[l3l[a)]. The fact 
that such a further variational improvement of the NRG 
ground state is possible, however, is of significance mainly 
as a matter of principle, not of practice: for the numer- 
ous situations where NRG works well (in particular, for 
A not too close to 1), we expect that such further vari- 
ational improvement of the NRG ground state will not 
noticeably affect any physical observables. 

Let us conclude with some comments about the pros 
and cons of NRG and VMPS. For quantum impurity 
models with a comparatively low degree of complexity, 
such as the single-lead Anderson and Kondo models, 
NRG works exceedingly well and for practical purposes 
nothing is to be gained from switching to the VMPS ap- 
proach. The latter is a potentially attractive alterna- 
tive to NRG only for two types of situations, namely (i) 
more complex quantum impurity models, and (ii) non- 
logarithmic discretization of the leads. We briefiy discuss 
these in turn. 

(i) For complex quantum impurity models, in partic- 



ular ones involving several leads, VMPS achieves a very 
significant reduction in memory cost, relative to NRG, 
for describing ground state properties via unfolding the 
Wilson chain. There are several caveats, though. Firstly, 
this reduction in memory cost applies only when only 
ground state properties are of interest. To obtain ex- 
cited state eigenspectra, the memory costs of NRG and 
VMPS are comparable. Secondly, unfolding is expected 
to work well only for models for which the subchains that 
are being unfolded are only weakly entangled, which will 
not be the case for all impurity models. For example, the 
two-channel model might be an example where unfolding 
works less well. In general, one needs to check the extent 
to which degrees of freedom on different subchains are 
entangled with each other, by calculating the mutual in- 
formation of two sites on different subchains. If this does 
not decrease rather rapidly with their separation from 
the impurity site, then unfolding will be a poor strategy. 
AppeaHngly, though, such a check can be done purely 
using NRG data, as illustrated in Section [III CI Thirdly, 
the fact that VMPS relies on variationally optimizing the 
ground state might cause convergence problems for mod- 
els which have degenerate ground states. Conceivably 
this problem can be reduced by systematically exploit- 
ing all relevant symmetries of the Hamiltonian, including 
non-Abehan symmetries,— >2^. However, if states in the 
local state space of a folded Wilson chain are related by 
a non-Abelian symmetry, then this symmetry would not 
be manifest in the unfolded representation. Thus, the 
two possible strategies for achieving significant memory 
reduction, namely unfolding and exploitation of symme- 
tries might not always be mutually compatible; which 
one is more favorable will depend on the details of the 
model, and is an interesting subject for further study. 

(ii) The VMPS approach offers clear advantages over 
NRG in situations where Wilson's logarithmic discretiza- 
tion of the conduction band cannot be applied. In the 
present paper, we found clear indications for this fact 
in the observation that the improvement of VMPS rel- 
ative to NRG becomes more significant as A is chosen 
closer to 1. More importantly, VMPS offers the pos- 
sibility, inaccessible to NRG, to improve the frequency 
resolution of spectral functions at high frequencies, by 
using a fiexible (non-logarithmic) discretization scheme 
which reduces the level spacing of effective lead states 
in the energy regimes where higher frequency resolution 
is desired. For such a discretization scheme Wilsonian 
energy scale separation is lost and NRG truncation can- 
not be applied. However, the ground state can still be 
found variationally, and spectral functions can be com- 
puted using projection operator techniques. In this fash- 
ion, it has recently been possible to calculate the spectral 
function for the Anderson model at large magnetic fields, 
B > Tk, and to resolve the split Kondo resonance with 
sufficient accuracy to reproduce the widths expected from 
perturbation theory in this regime. These developments, 
though, go beyond the scope of the present paper and 
will be published separately .^^i^° 



17 



Acknowledgments 

We gratefully acknowledge fruitful discussions Frithjof 
Anders, Theresa Hecht, Andreas Holzner, Ulrich Scholl- 
wock, Prank Verstraete and Gergely Zarand, and thank 
Frithjof Anders for constructive comments on the 
manuscript. This work was supported by the Spintronics 
RTN and the DFG (SFB 631, SFB-TR12, De-730/3-2). 
Financial support of the German Excellence Initiative via 
the Nanosystems Initiative Munich (NIM) is gratefully 
acknowledged. 



Appendix A: TECHNICAL DETAILS 

In this appendix, we collect some technical details on 
various manipulations involving matrix product states. 



1. Orthonormalization of B-matrices of unfolded 
Wilson chain 

To keep the notation simple, in this subsection we 
shall imagine the sites of the unfolded Wilson chain to 
be stretched along a line running from left to right, 
enumerated by an index k running from 1 for site Nl 
to K = 2{N + 1) for site iVf. Correspondingly, ma- 
trix product states will generically be written as j^*) = 

'E{,r'<} k^XrifcLi B^""^), with matrix elements ^i'^"'. 

It is convenient to ensure that every S-matrix in a 
matrix product state satisfies one of the following two 
orthonormality conditions: 

^^k-^lt^k.] ^ 1 ^ (Ala) 
^5k<.]sk<=]t = 1 . (Alb) 

In particular, if all i3-matrices satisfy either the first or 
the second of these conditions, the corresponding matrix 
product state is automatically normalized: 

(A2) 

This follows by iteratively applying Eq. (|Aip . To 
start the iteration, note that for matrices at the be- 
ginning or end of the chain, where one of the matrix 
indices is a dummy index with only a single value, 
Eqs. lAEl or dAibl imply B^^,l^^ B^^^^^ = S^'^ or 

J2aK bI^i"^ B^^J^^ = S^^,, respectively. In the NRG ap- 
proach, all ^-matrices naturally satisfy Eq. IjAlaP [cf. 
Eq. ©]. 

In the VMPS approach, it is convenient to ensure that 
during variational optimization sweeps, Eq. (jAlaP holds 



for all matrices to the left of the site, say k, currently be- 
ing updated, and Eq. (|Albp for all matrices to its right. 
Thus, after optimizing the set of matrices Sl'^'^l at site 
k, this set should be orthonormalized before moving on 
to the next site, such that it satisfies Eq. (|AlaP when 
sweeping from left to right (or Eq. (jAlbp when sweeping 
from right to left). This can be achieved using singular 
value decomposition [cf. Eq. lfT7|l1: Arrange the matrix 
elements of the set of matrices into a rectangular 

matrix carrying only two labels, with matrix elements 
Bpri = -Si'r)''' (or Bi,fj — sl^rj'^), by introducing a compos- 
ite index P — {ak,v) (or 77 = (cr^jTy)). Using singular 
value decomposition [Eq. (fTZl) ]. write this new matrix as 
B ~ USV'^ . Then rewrite the matrix product of two 
neighboring B-matrices as = ijkfelijkfc+i] 

(or SI'^fe-ilBl'^"! = ^['^fe-ili?!'^'^]), where the new matri- 
ces B are defined by 

fil^l = U^^ , B%'^'^ - (5V^i?["'=+il)^5 , (A3) 

(or Bf;^ = Vl, , = {B^^^-^^USU )• (A4) 

The property U'^U = 1 (or V^V = 1) ensures that the 
new set of matrices i3['^'=l at site k is orthonormal accord- 
ing to Eq. IjAlaP (or Eq. jAlbJ), as desired. Now proceed 
to the next site to the right (or left) and orthonormalize 
^['^fe+i] (or B^^k-i]-^ ^ijg same manner, etc. 

The above procedure can be used to orthonormaHze 
the matrices of a randomly generated matrix product 
state before starting VMPS sweeping. Likewise, during 
VMPS sweeping, each newly optimized matrix can be or- 
thonormalized in the above fashion before moving on to 
optimize the matrix of the next site. 



2. Refolding 

This subsection describes how to refold an unfolded 
matrix product state of the form 

l^"r,>u= k")(-B[""^l...S[""ilB[""T]...sKT])^^^ 

(A5) 

shown schematically by sites nl to n\ of Fig.[Hc). Its two 
indices will be treated as a composite index fi = {f.ri) 
below. The variational matrix product state |4'''^)u of 
Eq. lfT5|) discussed in the main text is a special case of 
Eq. (|A5p . with n = N and v — — 1 . The goal is to 
express Eq. (|A5|) as a Hnear combination, 

a 

(/? — {v, rj) is a composite index) of an orthonormal set 
of "refolded basis states" of the form of Eq. (|2T1) , 

1^-^), = J2 k")(B[^°lB[^il . . . B[^"l)i„ , (A7) 



18 



shown schematically by sites to n of Fig. [DJa) . To this 
end, we precede iteratively in n. We use singular value 
decomposition to iteratively merge, for every pair of sites 
n| and n\ of the unfolded chain, the matrices -B|,'^"^' and 

B^^,"^^ into a new set of matrices B^^r} for site n of the 
refolded chain, thereby trading the indices cr„|,(T„| and 
vrj of Fig.[2l[c) for the indices cr„ and a of Fig.[2l[a). This 
is to be done in such a way that the matrices B^"^"^ are 
orthonormal in the sense of Eq. ([6]) , and that for the first 
few sites their dimensions increase in a way analogous to 
those of the y4[°'"l matrices of NRG, starting from 1 x d 
at site n — 0. 

For the first iteration step, start with n = 0, make a 
singular value decomposition of the matrix product 



{x]')^x\' 



(AS) 



with /?' = {v' , 7]'), and use U'^ to define a new set of d ma- 
trices _B[°'''1 for site of the refolded chain, with matrix 
elements 3"^°} = U^^a'- The B^'^°^ have dimensions 1 x d 
(the dummy first index has just one value), and are by 
construction orthonormal in the sense of Eq. |[6]), since 
U°''U° = 1. Upon inserting Eq. ||A8| into Eq. ^El, the 
factor produces the first matrix factor B^'^°^ in the 
refolded state (|A7p . thus completing the first iteration 
step. For the second iteration step, contract the factors 
5°V°t with the factors B^'^^i'^ and ^['"itI in Eq. {ASI, fac- 
torize the result as U^S^V^^ and use to construct new 
matrices B^'^'-^ for site 1 of the refolded chain, etc. To be 
explicit, for general n, make a singular value decomposi- 
tion of the matrix product 
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Figure 14: (Color online) Graphical representation of the vari- 
ational equation used for cloning, Eqs. I|32p or l|A14p . drawn 
for the case /i —I, and assuming all matrix elements to be 
real. The upper part of the figure represents ^ -^^kjj c{G\G) c; 
it simplifies to B'^^^i [left hand side of Eq. (MM] 

upon real- 
izing that the parts in dashed boxes represent the left hand 
sides of Eqs. l|A12all and IIA12bll . and hence reduce to unity. 
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with composite indices a = {an, a'), an — {ani,an]), 
P = {v,ri) and f3' = {v',ri'). Then use lA^ to define a 
new set of orthonormal matrices i?!'^"! for site n of the 
refolded chain, with matrix elements b"^^,"} — ^sa- In this 
way one readily establishes that |^'"^)u can be written in 
the form of Eq. {All, with C;j^ (5"V"t)„^. 

The dimensions of the matrices B^""^ initially grow by 
a factor of d with each iteration step, until their dimen- 
sions are restricted by the number of possible values of 
the composite index (3, namely D'j^ , with D'^ given by 
Eq. IIS]). Thus, the bI"^"! have dimensions x DJ^+i, 
with L*^ = min(d", J, which leads to Eq. l[20l) . 



Let fc/i label the site to be updated and write the cloned 
state, which is assumed to be of the form (flSl) . as 



Here we introduced the shorthands 



(^^'■Ir = ( 



= ( B^'^'"'" 



(AlO) 



(xf^) = (^['^"Jl ...bKmJ") ^ (Alia) 



(Allb) 



3. Cloning 

This subsection gives some details of the cloning pro- 
cedure of Section IVI Al The goal is to solve the varia- 
tional Eq. l(32|) . which determines the B-matrices of the 
cloned state \G)c- As described in the main text, this can 
be done by sweeping back and forth along the unfolded 
Wilson chain, and updating one matrix at a time. 



for the products of matrices standing before or after the 
one of present interest in the unfolded Wilson chain, and 
the labels ki^i or kr^J,r label the sites just before or af- 
ter this site. Moreover, assume that all the S-matrices 
in Xi and Xr have been orthonormalized according to 
Eq. (|AlaP or (|Albp , respectively. (This can always be en- 
sured by suitably orthonormalizing each i3-matrix after 
updating it, see below.) These orthonormality relations 
immediately imply similar ones for the matrix products 
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just introduced: 

V (x^'') (xi;''V = s,,, . (Ai2b) 

Thus, the norm of \G)c can be written as 

e(G|G)e = -^5]s!,':rl^i3b^l, (A13) 

vv' 

where A/" is a normalization constant ensuring that the 
norm equals unity. 

Using Eq. (|A13p . the variational Eq. l(32|) readily re- 
duces to 



TTa 



where {cr'^} denotes the local indices of all sites except 
the index of site k\x, and we have assumed all A- 
and -B-matrices to be purely real (exploiting the time- 
reversal invariance of the present model) . This equation 
completely determines the new matrix i3[°''=f 1 in terms of 
the ^-matrices specifying the NRG input state |G)f and 
the i?-matrices of sites other than the present one, which 
had been kept fixed during this variational step. 

Having calculated sl'^fef'], it should be properly or- 
thonormalized, following the procedure of Eq. I|A3P or 
Eq. (|A4p . depending on whether we are sweeping from 
left to right or vice versa. In other words, use the singu- 
lar value decomposition US\>'^ of the new-found matrix 
^['^fcfi]^ to transfer a factor or US onto its right or 
left neighbor, respectively, and rescale this neighbor by 
an overall constant to ensure that the new state |G)c 
is still normalized to unity. This concludes the update 
of site fc/i. Now move on to its neighbor, etc., and thus 
sweep back and forth through the unfolded Wilson chain, 
until convergence is reached. 
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